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TECHNICAL NOTE 39^9 


A THEORETICAL STUDY OF THE EFFECT OF UPSTREAM 
TRANSPIRATION COOLING ON THE HEAT-TRANSFER 
AND SKIN-FRICTION CHARACTERISTICS OF A 
COMPRESSIBLE , LAMINAR BOUNDARY LAYER 
By Morris W. Rube sin and Mamoru Inouye 

SUMMARY 


An analysis is presented which, predicts the skin-friction and heat- 
transfer characteristics of a compressible , laminar boundary layer on a 
solid flat plate preceded by a porous section that is transpiration cooled. 
The analysis is restricted to a Prandtl number of unity and linear 
variation of viscosity with temperature . 

The local skin friction has been found to have a low value in the 
region of transpiration cooling and then to increase until it approaches 
the value for a completely nonporous surface asymptotically. The initial 
increase in local skin friction is rapid as half of the ultimate increase 
occurs in a distance beyond the porous region that is about 20 percent of 
the length of the porous region for all rates of injection. When the 
total coolant flow rate is kept constant and the porous length is varied, 
it is found that the average skin friction on a partially porous plate is 
slightly lower than that on a fully porous plate. 

The local heat transfer behaves in a manner si m ilar to that of the 
local skin friction. It is found, in an example, that the temperature 
at the end of a partially porous plate could be maintained at about the 
same temperature as a fully porous plate by doubling the total rate of 
coolant flow. 


INTRODUCTION 


For flight at high speeds, aerodynamic heating often requires the 
cooling of aircraft in order to maintain tolerable surface temperatures . 

Of the various cooling techniques available, transpiration cooling systems 
are us ually effective for this application, as is shown in reference 1. 
This results because the geometry of the porous surface provides for 
excellent heat exchange between the coolant and the surface , and the 
bo undar y layer on the surface is altered so as to reduce significantly 
the skin friction and the heat transfer to the surface. 
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The attractiveness of transpiration cooling of large surfaces is 
reduced "by the introduction of structural problems. It is difficult to 
manu facture large porous surfaces and to support them in use because of 
their inherent weakness. To take advantage of transpiration cooling and 
also to alleviate the structural problems, the use of partially porous 
surfaces offers possibilities. In this scheme, the most critical regions 
from an aerodynamic heating standpoint could be transpiration cooled, and 
the downstream regions protected by the film of coolant that is introduced 
into the boundary layer. 

It is the purpose of this investigation, therefore, to examine 
theoretically the magnitude of protection offered the downstream nonporcrus 
regions by the film cooling process. The present analysis is restricted 
to consideration of a compressible, laminar boundary layer on a semi- 
infinite flat plate with air as the coolant. The flat plate is divided 
into two regions: an upstream region of finite length which is porous 

and transpiration cooled, and a downstream region which is nonporous and 
is protected by the upstream cooling process . 


SYMBOLS 


a n (x), a 


a B 


C 


Cf 

Cp 

Do 


f (0) = F 


k 

L 

Rr 


dimensionless coefficients in assumed velocity profile 
polynomial (A19) 

(The subscript is dropped for n = 1.) 

constant value of a for no transpiration or Blasius solution 

constant of proportionality between absolute viscosity and 
absolute temperature defined in equation (A 6) 


local skin-friction coefficient. 


t s 


1 2 
g Poo^oo 


specific heat at constant pressure 


R 


constant, 


So 




dimensionless stream function used in reference 2 such that 
f'(tl) - 


dimensionless number which is proportional to the mass flow 
of transpiration cooling 
(See eq. 5.} 

thermal conductivity of fluid 
length of plate 

C'nLL 

Prandtl number, 
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P 

CL 

R 

^x: 

R S> 

R 8 q > r L 
s(x) 

0? 

u 

V 

w 

X 

X' 

y 

y* 

r 

8 

8+ 

T1 

^8 

9 

P 

v 


pressure 

heat-transfer rate per unit area 
gas constant 

Reynolds number, based on distance 


' CPco 


Reynolds numbers based on 8, x Q , 5 0 , L 

slope of velocity profile at surface, 
absolute temperature 
velocity component parallel to surface 
velocity component normal to surface 
total coolant flow rate 

transformed coordinate along surface defined in equation (A9) 
physical coordinate along surface 

transformed coordinate normal to surface defined in 
equation (A9 ) 

physical coordinate normal to surface 
gamma function 

boundary-layer thickness in transformed coordinates 

8 Rs 

dimensionless boundary-layer thickness, -s— = = — 

5 ° R S 0 

similarity parameter used in reference 2, 

value of T) corresponding to the edge of the boundary layer 

dimensionless temperature function defined in equation (A59 ) > 
T-5 s 

^o-T s 

absolute viscosity 
kinematic viscosity 
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P 
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9(e) 


c 

cl 

s 

t 

00 

O 


( ) 


dimensionless length along surface, — L = 


Rv 


L o ^x 0 


density 

dimensionless slope of velocity profile at surface, 
a 80 / x 


shear stress at surface 

function defined in equations (A83) 


Subscripts 


transpiration cooled surface 

coolant 

surface 

stagnation conditions 
free stream 

conditions at end of porous region 


Superscripts 

quantities related to original physical coordinates 
average 


ANALYSIS 

Conditions and Assumptions 


The analysis is performed for a compressible, laminar 
boundary layer on a nonporous flat plate behind a porous region -which is 
transpiration cooled. The geometry is shown in figure 1. The analysis 
can employ any erf existing solutions (such as refs. 2, 3, and 4) for the 
transpiration cooled region, and use their results at x* = x 0 ’ to begin 
the present analysis. For convenience, the solution (refs. 2 and 3) 
containing the following conditions is employed: 

1. The coolant is the same as the boundary-layer fluid. 
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2. The velocity of the coolant normal to the surface is inversely 

proportional to the square root of the distance from the leading 
edge . 

3 . The porous surface temperature is constant. 

In addition, the following assumptions concerning the boundary-layer 
fluid are made: 

1. Brand tl number is unity. 

2. Specific heat at constant pressure is constant. 

3- The fluid behaves as a perfect gas; p = pRT. 

1+. The absolute viscosity varies linearly with the absolute 
temperature . 


Summary of Analysis 


The mathematical details of the analysis are included in Appendix A, 
and only a summary is presented here. 

Transformation of basic equations .- The analysis is begun with the 
basic equations of continuity, momentum, and energy for a compressible, 
laminar boundary layer with no pressure gradient. By use of the 
Stewartson transformation (ref. 5 ) in which the physical coordinate normal 
to the surface is modified to take into account the compressibility, the 
basic equations are reduced to the equations for an equivalent, incom- 
pressible flow. This transformation not only simplifies the equations , 
but together with the previous assumptions concerning the fluid, permits 
solution of the momentum equation independently of the energy equation. 

Solution of momentum equation for skin friction .- To begin the 
present solution, it is necessary to start with velocity distributions 
that match those of the boundary layer at the end of the porous region 
where transpiration cooling ceases . Because these velocity distributions 
change from those characteristic of transpiration cooling to those char- 
acteristic of flow along a solid surface far downstream along the plate, 
similarity of profiles cannot be expected, and an exact solution is not 
simple to obtain. Therefore, the momentum integral method of solution 
was chosen. 

A convenient and fairly accurate velocity distribution can be found 
using the following procedure. The boundary-layer velocity profile for 
use in the momentum integral equation is approximated by a seventh-degree 
polynomial in terms of y/8. 
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where y is the transformed distance normal to the surface and 8 is 
the transformed boundary- layer thickness . The coefficient of the first 
term, a, is a measure of the shear stress or the slope of the velocity 
profile at the surface. The seven unknown coefficients of the polynomial 
ft.nfl the boundary-layer thickness make a total of eight unknowns to be 
determined. In addition to the momentum integral equation, seven other 
equations are necessary in order to solve the problem. These equations 
are obtained by imposing seven boundary conditions on the velocity profile. 
Four of the conditions are imposed at the outer edge of the boundary layer 
and three at the surface. To insure that the local velocity in the 
boundary layer approaches the free-stream velocity smoothly, it is required 
that for y = 5 




U = -On 


du _ ^u n 

by dy 2 Sy 3 

• 

To insure that the velocity profile has the correct slope and curvature 
at the surface, the velocity distribution and the boundary conditions 
u = v = 0 are inserted into the basic momentum equation and its deriva- 
tives with respect to y. This yields for y = 0 (see eqs. (Al6), (A2^), 
and (A25 ) ) 


e^u _ ^u _ o 
by 2 dy 3 

du 6 /chA _ „ d 4 u 
by ctx \^y/ 00 Sy 4 ( 


(3) 


These boundary conditions permit the evaluation of the unknown coefficients 
of the polynomial in terms of the first coefficient a, and lead to an 
ordinary differential equation relating a, 5, and x. This equation and 
the momentum integral equation are solved simultaneously along the plate. 

Before the differential equations can be solved, initial conditions 
for these differential equations must be determined. The initial condi- 
tions for the present solution are obtained by matching the seventh-degree 
polynomial velocity profile to the exact velocity distribution at the end 
of transpiration cooling. This is accomplished by equating the slopes at 
the wall and the boundary-layer thicknesses of the two velocity distribu- 
tions. A problem arises in the definition of the boundary- layer thickness. 
In the exact solution, the local velocity in the boundary layer approaches 
the free-stream velocity asymptotically, whereas with a polynomial solu- 
tion the local velocity equals the free-stream velocity at a finite 
distance from the surface. In this analysis, the boundary-layer thickness 
for the exact velocity distribution is chosen as the point where 
u/uoo = 0.9976, and the reason for this choice is in Appendix A. 
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With these initial conditions, the system of two s imult aneous 
ordinary differential equations was integrated numerically yielding 5 
and a as functions of x. The local skin-friction coefficient was then 
calcinated directly from these results. 

Solution of energy equation .- Normally, the surface temperature of 
the nonporous region will vary and differ from the constant surface tem- 
perature of the transpiration cooled region. Hence, the solution of the 
energy equation must account for surface temperature variations. A con- 
venient method of analysis is to separate the solution of the energy 
equation into two parts. First, a solution of the complete inhomogeneous 
energy equation is found for the condition that the porous and nonporous 
regions are at the same constant temperature. Then solutions of the 
homogeneous portion of the energy equation are found and added to satisfy 
the surface temperature boundary conditions. This addition of solutions 
is permissible because the energy equation is linear in temperature. 

The inhomogeneous energy equation can be expressed in terms of total 
energy such that it has the same form as the momentum equation (since 
Pr = l). Then, for the ease of constant surface temperature, the boundary- 
layer tenperature distribution can be expressed in terms of the velocity 
distribution, which has been solved previously. 

Initially, the homogeneous portion of the energy equation is solved 
for the condition of a single surface temperature discontinuity behind 
the transpiration cooled region. In order to obtain a solution, the 
boundary-layer velocity profile Is replaced by a linear profile having 
the correct slope at the surface as determined previously in the solution 
of the momentum equation. Then the temperature distribution In the 
boundary layer can be solved for In terms of a single similarity 
parameter, £. 

The general solution of the energy equation, which takes into account 
arbitrary surface temperature variations, is obtained by adding the solu- 
tion for the inhomogeneous energy equation to as many solutions of the 
homogeneous equation as required to satisfy the surface temperature 
boundary condition. For an arbitrary surface temperature distribution, 
the summation of the solutions of the homogeneous energy equation can be 
expressed as an integral. Hence, the general solution of the energy 
equation yields the temperature distribution in the boundary layer for 
an arbitrary surface temperature distribution. 

The heat-transfer rate at the surface is obtained by differentiating 
the temperature distribution with respect to y*. If the surface tem- 
perature is given, the heat transfer can be determined. For the case of 
a constant surface tenperature, the Reynolds analogy is applicable, and 
the local heat transfer is directly related to the local skin friction. 

Thus far, it has been assumed that the surface tenperature distribu- 
tion is given as a boundary condition. In many cases, however, the heat 
transfer at the wall is given as a boundary condition, and the surface 
tenperature distribution is desired. This problem is solved by using an 
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integral equation based on the equation for the heat-transfer rate. In 
particular, the surface temperature distribution was solved numerically 
for the condition that the nonporous surface is insulated. 


Results of Analysis 


The results of the analysis are summarized in the following equations 
The local skin friction is given by 


ff 

2 



a 

B + 





where 


t 


X 

x O 


Note that both a(l) and Rg o are dependent on the transpiration coolant 

rate in the porous region and the Reynolds number of the porous region. 
For example, when the local blowing rate is expressed by f(0) = F where 


f(0) = 


-2 


Ps^s 


• 7 % = - 


Ps v s rrr* 


(5) 


then 


R &o 




( 6 ) 


c(l) = a(F) 


(7) 


Note that PsV s = 2p s vs because Vg varies inversely as For 

convenience D 0 (F) and cr(F) are tabulated as functions of F = f(0) in 
table I. For § > 1 or x > Xq a plot of cf(|) for each f(0) is given 
in figure 2. 

If the surface temperature distribution is given, the local heat- 
transfer rate at the wall for x > Xq is found from 
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is(s) = -Z^-RT) 

R 6 0 


(V T c> v/a( |) -0.538Do 37 a y 


I 


dJ B (Si) 
as a 


46 , 


[«P( 6 )-<P( 6 1 )] 


i/3 


( 8 ) 


where 



(9) 


The function cp(|) is presented in figure 3 for various f(O). 

If the heat-transfer rate at the surface is given, the surface 
temperature distribution can he determined from 


* S (§)-T C 


0,512 r 6 ^w -' + cvtoMU 

Do 1 ' 3 '! [^D-tpU ,.)] 2 ' 3 


( 10 ) 


RESULTS AHD DISCUSSION 

Comparison of Polynomial and Exact Velocity Distribution 


As explained in the analysis, at the start of the nonporous region, 
the seventh-degree velocity profile must match the exact transpiration 
cooling profile as given in reference 2. To do this, the polynomial 
velocity was made to have the same . slope at the surface and to have the 
same boundary-layer thickness as the exact velocity profile. A comparison 
of the two entire profiles is shown in figure k for three tran s piration 
cooling rates where the absolute magnitude of f ( 0 ) is proportional to 
the amount of transpiration cooling. The agreement between the reppe&tive 
profiles throughout the boundary layer -is satisfactory. Because the 
Prandtl number is assumed to be unity, (agreement of the velocity profiles 
implies an agreement of the corresponding temperature distributions. 

Skin Friction • 


The effect of upstream surface transpiration cooling on the local 
skin friction is shown An figure 5« The abscis&a Is' the ratio of the 
distance from the leading edge to the length of the porous region. The 
ordinate is the local skin-friction coefficient with upstream .transpiration 
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cooling divided by the local skin friction that would exist with no 
transpiration cooling. This figure was obtained from figure 2. For 
each transpiration cooling rate, the reduction of skin friction diminishes 
rapidly with downstream distance and approaches the no transpiration solu- 
tion asymptotically. For example, continuous transpiration cooling of 
f(0) = -0.50 halves the skin friction. With a finite transpiration-cooled 
region, however, the skin friction rises to within 13 percent of the no 
transpiration value at 5=2 and to within 5 percent at | =6. Similar 
effects are shown for the other transpiration coolant rates. 


•Heat Transfer 


The effect of upstream surface transpiration cooling on the heat- 
transfer characteristics of the nonporous region was evaluated for two 
cases: the local heat transfer at the surface was determined if the 

surface temperature distribution was prescribed, and the surface tempera- 
ture distribution was determined if the local heat transfer at the surface 
was prescribed. 

The first case occurs when the surface temperature is constant over 
both the porous and nonporous regions. By Reynolds analogy, the effect 
of upstream surface transpiration cooling on the heat-transfer rate is 
the same as the effect on the skin friction shown in figure 5- For 
example, at 5=2 with f(o) = -0.50, the internal heat-transfer rate 
necessary to maintain a constant surface temperature is reduced 13 percent 
from the no transpiration value. Similar trends occur for the other 
transpiration coolant rates. 

An example of the second case occurs when the nonporous region is 
insulated, or qg(5) = 0. The resulting surface temperature distribution 
is shown as the ordinate in figure 6. The ordinate is written in dimen- 
sionless form employing the temperature of the cooled upstream region and 
the total temperature. It is observed for an coolant rates that the 
surface temperature of the nonporous region rises rapidly at first and 
ultimately approaches the recovery temperature asymptotically. The 
reduction in temperature due to the effect of injecting transpiration 
coolant into the boundary layer is not very large, however, as is noted 
when the rest of the curves are compared with the f ( 0 ) = 0 curve, which 
represents the temperature distribution caused by cooling the upstream 
section to T c by some internal cooling system. 


Practical Implications of Results 


Figures 5 and 6 show how the skin friction and heat transfer behave 
on the nonporous section when air is transpired in the porous section. 
The parameter “of these curves, f(0), however, does not provide a satis- 
factory criterion for a practical evaluation of the usefulness of the 
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cooling system under analysis . A more useful criterion is the total 
coolant flow rate. This permits comparison of various cooling systems on 
a corresponding weight basis. 

For instance, it is interesting to know how the average skin friction 
on a plate changes as the proportion of the lengths of the porous and 
nonporous regions is varied, subject to the restriction that the total 
coolant flow rate be constant. Thus, as the porous region is reduced in 
length (actually area), the rate of coolant flow per unit area is 
increased. It is quite easy to calculate these effects using the results 
of figure 5* Thus, in figure 7 > the ratio of the average skin friction 
on the plate to the corresponding average skin friction on an entirely 
solid plate is plotted as a function of the fractional length of the plate 
that is porous. The parameter of each curve is a dimensionless group con- 
taining a term representing the total coolant flow rate, w. Starting from 
the right side of the figure, it is noted that for a completely porous 
plate (xq/L =1) there is a considerable reduction in skin friction for 
each coolant rate. As the fraction of porous surface is reduced, x 0 /L 
becoming s malle r, it is rather surprising that the reduction in average 
skin friction is even greater for each of the total coolant rateB. Appar- 
ently, the greater reduction in the local skin-friction coefficient over 
the porous section of the plate due to the higher local coolant rate more 
than compensates for the rise in local skin friction experienced over the 
nonporous section as shown in figure 5* The line, at the left, represents 
an estimate of the lower bound of these results, because at a fixed total 
coolant rate the local coolant rate at points to the left of this dashed 
line becomes sufficiently large to separate the laminar boundary layer. 

The curves are drawn as far to the left as the information in figure 5 
allows. It appears, therefore, that on the basis of average skin friction, 
no penalty results from restricting the transpiration cooled portion to a 
fraction of the total length of the plate when a given total amount of 
coolant Is used. 

In figures 8(a) and (b), there are shown the temperature distributions 
that result from using a fixed quantity of coolant flow with various 
lengths of porous regions. The ordinate represents the surface tempera- 
ture expressed in terms of the recovery temperature, which for Pr = 1 
is the total temperature, and the initial coolant temperature. It can be 
seen in figure 8(a) that for the amount of coolant used, a fully porous 
plate would result In a dimensionless surface temperature of about 0.66. 

The temperature distributions resulting from the plate being 0.5, 0.25, 
and 0.1 porous are indicated also. For these cases, the temperature of 
the porous regions is reduced markedly because of the high local rate of 
coolant flow. The reductions in the porous surface temperature, however, 
are not sufficient to lower the surface temperature over most of the rest 
of the surface down to the value of the fully porous case. It is signi- 
ficant that even at the end of the plate, in these examples, the surface 
temperature is well below the recovery temperature. As an example of 
these temperatures consider an aircraft flying at a Mach number of 5 
under steady-state conditions, where the recovery temperature is about 
1770°F. If cooled as in figure 8(a) with an entirely porous surface and 
a laminar boundary layer, the surface temperature would be about 1220° F 
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for an initial coolant temperature of 150° F. With partial cooling, the 
maximum temperature at the trailing edge would be about l460° F, or still 
about 300° F below the recovery temperature . 

The effect of doubling the coolant flow rate is shown in figure 8(b). 
For the fully porous case the dimensionless surface temperature becomes 
about 0.4. The only partially porous case shown is where the plate is 
0.25 porous. Again the dimensionless temperature is low over the porous 
region (about 0.07) and rises rapidly in the nonporous region reaching a 
value of 0.66, which happens to be the value obtained on the fully porous 
plate with half the coolant flow. Thus, to maintain a prescribed maximum 
surface temperature , the partially porous surfaces will require more total 
coolant than fully porous surfaces. 


CONCLUDING REMARKS 


An analysis has been made that determines the effect of upstream 
surface transpiration cooling on the skin-friction and heat-transfer 
characteristics of a compressible, laminar boundary layer on a flat plate. 

The skin friction has been found to have a low value in the region 
of transpiration cooling and then to increase until it approaches asymp- 
totically the value for a completely nonporous surface. The initial 
increase in skin friction is rapid as half of this ultimate increase 
occurs in a distance beyond the porous region that is only about 20 percent 
of the length of the porous region for all rates of injection. When the 
total coolant flow rate is kept constant, however, it is found that the 
average skin friction on a partially porous plate is even lower than that 
on a fully porous plate. This occurs because the local transpiration rate 
per unit area increases for the shorter porous regions so that additional 
reductions in skin friction in the porous region compensates for the rapid 
increase in the skin friction over the nonporous region. 

Heat-transfer solutions for the nonporous region were obtained for 
two types of problems: the local heat-transfer rate can be found if the 

surface temperature distribution is given, and the surface temperature 
distribution can be found if the local heat-transfer rate is given. A 
detailed example was presented - the solution to the problem of determin- 
ing the temperature distribution along the nonporous surface so that there 
is zero heat transfer at the surface. This corresponds to a thin surface 
at equilibrium behind a transpirationr cooled region. From this example, 
it was found that the temperature rises from the value of the transpira- 
tion cooled surface and approaches the recovery temperature at large 
distances along the plate. The initial rise in temperature is quite 
rapid as half of the ultimate rise in temperature occurs in the region 
between 1.3 and 2.0 times the porous length, for the range of transpira- 
tion coolant rates considered. For constant coolant flow rates, with 
varying proportions of the porous region length to the total length, it 
is found that the additional cooling in the porous region due to the 
higher local coolant flew does not compensate for the rapid rise in 
temperature over the nonporous region so that the surface temperatures 
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at the end of the plate are larger than would exist on a f ully porous 
plate. For the examples shown, the temperature at the end of the par- 
tially porous plate could "be maintained at about the same temperature as 
a fully porous plate by doubling the total rate of coolant flow. In 
practical application, the increased flow requirements of the partially 
porous system may be offset somewhat when considerations are mad e of the 
weight of the porous surface supporting structure and of thermal radiation 
from the surface on transient effects which normally cause the rear 
regions on a surface to have less temperature problems. 


Ames Aeronautical Laboratory 

National Advisory Committee 
Moffett Field, Calif., 


for Aeronautics 
March 7, 1957 



NACA TN 3969 


14 


APPENDIX A 

MATHEMATICAL DETAILS OF TEDE ANALYSIS 
Transformation of Basic Equations 


Basic equations, assungtions , and Boundary conditions .- The "basic 
equations for a compressible, laminar boundary layer with no pressure 
gradient are: 

Continuity equation 


a 

ax' 


(p , u , )+ 


a 

aF 


(p*v») 


0 


Momentum equation 


p’u* 



p’v* 



auA 
Sy v 


Energy equation 


aT* ar* 

P ,u,c pf7 + P^ ,c pf 7 



(Al) 


(A2) 


(A3) 


A number of assumptions about the boundary-layer fluid are made to 
simplify the analysis: 

1. Pr = . 1 (Alt) 

2. Cp = constant 

3. The fluid is assumed to behave as a perfect gas, which has an 
equation of state. 


p* = p*RT* 


(A5) 


b. 


al = cSl 

Poo IS 


(A 6) 


As a result of the first two assumptions, the energy equation (A3) 
becomes 


p'u> 


' 


ar 1 

ay* 



(AT) 
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The boundary conditions for the nonporous region, x* > x 0 ’, are 


y t =o: u* = 0, v* = 0, T* = T s »(x») I 


y* -^00: T»-> Tx, | 


(A8) 


Stewartson transformation of variables .- The independent variables 
(x 1 , y*) are replaced by a new set of variables (x, y) in which the y* 
coordinate is modified to take into account the compressibility as done 
by Stewartson (ref. 5). 


x 


y 



(A9) 


The formulas relating the derivatives with respect to the physical 
and transformed independent variables are 


f d \ 

WVx' Poo W/x 

(hi, ’ (k)/(hl,(kl 


(A10) 


The dependent variables are then related to a stream function ijr, 
which satisfies the continuity equation and is related to the velocity 
components as follows : 


P*u’ = fbo 


in 1 

3y* 


p ' T ' = h 


(All) 


Substitution of equations (AlO) and (All) with the assumptions (A5) 
and (A6) into equations (A2) and (A7) yields 
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(A12) 


'tXrMtXf) 


= Cv„ 



(A13) 


To put these equations in the familiar form of the incompressible 
equations, let u and v be defined by 


u 


M 1 

Sy j 


v 


& 

J 


(Alh) 


It is noted from equations (A10) and (All) that u = u’, or that the 
transformed velocity and the physical velocity components in the x or x* 
directions are identical. The preceding definition for u and v results 
in a continuity equation identical to that for an incompressible flow. 


du 

8 x by 



(A15) 


Substitution of equations (A14-) into (A12) and (A13) results In 
momentum and energy equations Identical to those for an incompressible 
flow. 


u 


5u 

u +v 
dx 


&T St 


Su 

¥ = 


CVoo 


CVo 


S 2 ! 


SSu 

(ai 6 ) 

+ i(i)] 

(A17) 


The prints on the T l s are dropped for simplicity. Hence, by use of the 
Stewartson transformation, the compressible boundary-layer equations have 
been reduced to equivalent incompressible equations. 

In terms of the new coordinates, the boundary conditions (A 8 ) for 
the nonporous region, x > x OJ become 
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y = 0: u = 0, v = 0, T = T B (x) 

y-> oo; u-> u», T-> 


(A18) 


Solution of Momentum Equation for Skin Friction 


Seventh-degree polynomial velocity profile and evaluation of 
coefficients The houndary-layer velocity profile is approx ima ted hy a 
seventh-degree polynomial in terms of y/6, where & is the transformed 
houndary-layer thickness. 


u 

Uoo 


n=r 

= ^ an(x)(p n 
n=i 


(A19) 


The mth derivative with respect to y is 

n=T 

(A2o) 


The slope at the surface is 

(A21) 


where a is the coefficient of the first power term and is a measure of 
the shear stress at the surface. 

The introduction of the seventh-degree polynomial velocity profile 
adds seven unkn own coefficients to the problem, in addition to the 
houndary-layer thickness, for a total of eight unknowns . With the momen- 
tum integral equation as one equation, seven other equations are required 
in order to obtain a solution. These equations are obtained hy imposing 
seven conditions on the seventh-degree polynomial velocity profile. Four 
of the conditions are boundary conditions imposed on the velocity and its 
derivatives at the outer edge of the boundary layer. To insure that the 
local velocity in the boundary layer approaches the free-stream velocity 
smoothly, it is required that for y = &, 
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U = Uoo 

8u _ n 

8y Sy 2 


(A22) 


Substituting these conditions into equations (A19) and (A20) yields 


a+as+as+a^+as+ae+ay = 1 
a+2a2+3a 3 +4a4+5a5+6a6+7a7 = 0 
s>2+3a3+6a4+10aB+15a6+21 a 7 = 0 
a3+4a4+10a5+20ae+35a7 = 0 


(A23) 


Before the remaining conditions are imposed, the first and second 
derivatives with respect to y of the momentum equation (Al6) are 
required. The first derivative is 


^u _ 

SxSy Sy 2 


CVoo 


8 2 U 

Sy 3 


(A2^) 


and the second derivative is 


^u 8 2 u + u 8^ 
8y 8x8y dxdy 2 


8u a^u 

8x Sy 2 


+ v 



Cv Q 


a^ 

ay 4 


(A 25 ) 


The three remaining conditions are imposed at the surface. To insure that 
the velocity profile has the correct slope and curvature at the surface, 
the seventh-degree polynomial with the boundary conditions, u = v = 0, is 
substituted into the momentum equation (Al6) and its derivatives (A2^) 
and (A25). There results 


a£ — £13 = 0 


(A26) 


and an ordinary differential equation 
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Equations (A23) and (A26) constitute a system of six algebraic 
equations with seven unknowns. These equations are solved simultaneously 
to obtain a] 1 the unknowns in terms of the coefficient of the first 
term, a . 


&2 = 0 
&3 = 0 

a 4 = 35 - 20 a 
a 5 = 45a- 84 
as = 70-36a 
a T = 10 a -20 




(A28) 


Substituting a 4 from equation (A28) into equation (A27) and 
carrying out the differentiation yields 


IS 2 ( 5 S' : a i) - (35-20a)24 (A29) 


This equation, in part, determines how the slope of the velocity profile 
at the surface varies with the distance along the surface. 

Momentum integral equation .- The momentum integral equation is 
obtained by integrating the momentum equation (Al 6 ) from y = 0 to y = 5. 

CVqq 
Uco 2 



The assumed velocity profile (A19) with the coefficients expressed in 
terms of a from equations (A28) and the slope at the surface given by 
equation (A21) are substituted into equation (A30). The integration is 
carried out to yield 


- 4 E 5 ) ] 


(A31) 


where 


Aj. = 0.09518 
A 2 = 0.04371 

A3 = 0.01832 




(A32) 
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Reduction of problem to solution of two simultaneous ordinary 
differential equations Equations (A29) and (A31) form a system of two 
simultaneous ordinary differential equations in a, 5, and x. These 
equations can he simplified by expressing S and x in terms of Reynolds 
numbers • . 


R 5 Cv, 


&Uco 
'00 


R-sr = 


_ XUco 


CVo« 


(A33) 


Substituting these Reynolds numbers into equations (A29) and (A31) 
yields 

^ - a - 2M35-20a) <«>0 


^ = d^^ R s( A 2- +A 2a'-A3a 2 )] 


(A35) 


The preceding equations can be manipulated to separate the derivatives 
and yield 


dRs _ a x (a) 

&R X R& 


(A36) 


da _ cxg(a ) 
&R X " R§ 2 


(A37) 


where 


01 (a) = 


/ \ 120 
02 (a) = -g- 


7Ai+(7A 2 ' 


“7A 2 +( 4a 2 +14a 3 ) a +^^Q - 8 A. 3 ^a 2 

Ai+2A 2 a-3A3a 2 u 

r k A 1 )a-(4A 2 +7A 3 )a 2 +^—^ +4A: 


las 


A x +2A 2 a- 3 A 3 a c 


The Reynolds numbers Rg and R x can be normalized by substituting 
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Jk 

R X„ 


x 

xT 


S+»Jh = 8. V 

*5„ 8 c / 


D n = 




Equations (A 36 ) and (A3T) then become 


dS+ a,i(a) 


d| D 0 5 + 

da _ 02 (a) 

d| " D 0 (&+) S 


(A39) 


(a4o) 

(a4i) 


The initial conditions at this stage are 


| = 1 , 5 + = 1 


(a42) 


Before the equations can be solved, it is necessary to know the initial 
conditions for a and D 0 * 

Determination of initial conditions for solution of momentum 
equation .- Before a solution can be obtained for equations (A4o) and (a 41 ) , 
it is necessary to obtain the initial conditions for a and D 0 by match- 
ing the seventh-degree polynomial, both in thickness and slope at the 
surface, to the exact profile at the end of the transpiration-cooled 
region. Since the seventh-degree polynomial profile assumes a finite 
boundary-layer thickness, whereas in the exact profile the velocity 
approaches the free-stream velocity asymptotically, a problem arises in 
defining the outer edge of the boundary layer. In this analysis the no 
transpiration or Blasius case is used to make thiB definition. 

For the Blasius case, the shear stress at the surface, or the skin- 
friction coefficient, varies as R x “ 1/2 , and the boundary-layer thickness 
varies as R x 1/2 . The product of the two quantities is, therefore. Inde- 
pendent of Rx- For the present solution, the product of the shear stress 
at the surface, which Is proportional to the slope at the surface as given 
in equation (A21), and the boundary-layer thickness is a f unction of a 
alone. For a to be independent of R x , it follows from equation (A37) 
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that 02 (a) =0 or ag= 1.871. This value of a B for the seventh-degree 
polyno mial with no transpiration cooling must now he related to the exact 
solution. 

The following relations for the exact solution are obtained from 
reference 2. 


% = ^ >J ^x 

Cf _ 1 f f, (0) 

~2 ~ k 

15) 


(A43) 

(a44) 

(a45) 


The function f (-q ) is a dimensionless stream function, and rj is the 
similarity variable with r\§ denoting its value at the outer edge of the 
boundary layer. 


The local skin-friction coefficient is defined a6 


2 fto^xj 


(a46) 


where 


T 


I 

s 



(A4T) 


For the present analysis, the local skin-friction coefficient is 
obtained by substituting equations (A5), (A6), (AlO), (A2l), and (A47) 
into equation (A46). There results 


2f 

2 



(A48) 


The boundary-layer thickness and the local skin-friction coefficient, 
which is determined by the slope of the velocity profile at the surface, 
for the seventh-degree polynomial solution and the exact transpiration 
cooling solution are now equated. From equations (a43), (A44), and (a48) 
there results 


_ _ 2a 


(A49) 
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With no transpiration cooling, ag = I. 871 , f"(0) = 1.3282, and 
7]§ = 2.SL8. From table II of reference 2 and the relation (A^5 ) , the 
outer edge of the boundary layer is defined as the point where the 
velocity ratio is 

— = | f *Cns) “ 0.9976 (A50) 

This definition is used to define the boundary-layer thickness when 
transpiration cooling occurs. 

The initial conditions for the solution of equations (A4o) and (A4l) 
are thus determined as follows: A value of f(0) is selected using 

equation (5), the absolute magnitude of f(0) being proportional to the 
amount of transpiration cooling. With the outer edge of the boundary 
layer defined by equation (A50), the value of T]g is determined frcan 
table II, reference 2. Then a is calculated from equation (A^9), and 
D 0 is calculated, from equation (A^3 ) . With these initial conditions, 
the system of two simultaneous ordinary differential equations (a4o) 
and (A4l) was integrated numerically on an IBM 650 computing machine to 
obtain a and 8+ as functions of g for £ = 1 to £ = 10. Five differ- 
ent transpiration cooling rates were selected, including the Blasius case. 
The results are expressed in terms of the local skin-friction coefficient 
as given by equation (A48). 


Solution of Energy Equation 


Separation of solutions . - Because the surface temperature of the 
nonporous region will normally differ from the constant surface tempera- 
ture of the transpiration-cooled region, it is necessary in this analysis 
to obtain a solution of the energy equation (A17) that accounts for sur- 
face temperature variations . The method of solution followed will be to 
obtain a solution for the complete inhomogeneous energy equation when the 
porous and nonporous regions are at the same temperature . Then solutions 
of the homogeneous portion of the energy equation are added to satisfy 
the surface temperature boundary conditions . This addition of solutions 
is permissible because the energy equation is linear in temperature. 

Solution of inhomogeneous energy equation .- The inhomogeneous energy 
equation (A17 ) is solved for the case of a constant surface temperature , 
T c . If the momentum equation (Al6) is multiplied by u/cp and added to 
the energy equation (A17), there results 

u s ( T+ £f) +v i ( T+ i|) - ( T+ s|)] (A51> 

Because this equation has the same form as the momentum equation, if u 
is a solution of the momentum equation, then 
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T+~ = CiU+C 2 
2c p 


(A52) 


must be a solution of the energy equation. When the boundary conditions. 


u = 0, T = T c 
u = Uoo, T = %o 


(A53) 


are imposed, the unknown constants. 


Cl 

0 Z = T c 


(A?h) 


are determined. Hence, an exact solution of" the energy equation for the 
case of a constant surface temperature is 

T I - <V*o>& + Vsi| («5) 

Solution of homogeneous energy equation for single discontinuity in 
surface temperature .- The homogeneous energy equation is solved for a 
single discontinuity in surface temperature. The boundary conditions on 
the surface temperature are for 


x < xj , T = ^ 


x > Xj , T = T s 


(A 56) 


where xj > x 0 is the location of a surface temperature discontinuity. 

In order to obtain a solution of the homogeneous portion of the 
energy equation, which is equation (A17) with the viscous dissipation 
term neglected, the boundary-layer velocity profile is assumed linear 
with the correct slope at the surface, as determined previously in the 
solution of the momentum equation. The x component of velocity is 
assumed to be 


u = s(x)y 



(A57) 
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The y component of velocity is obtained by integrating the continuity 
equation (A15) with respect to y. 


v = 




(A58) 


The temperature distribution in the boundary layer is assumed to be 
a function of a single parameter £, where 


T-T5 

^50 _ f S 


= see) 


(A59) 


and 


e = X(x)Y(x, xj)y (A60) 


where X and Y are two arbitrary functions of x and x^. From 
equations (A 56 ), (A5 9), and (A60 ) , the boundary conditions on 0 are 


X < XjJ 

all y, 

0=1 



all x: 

y-> 

0=1 

1 

(a 6 i) 

X > Xj. 

y = 0 , 

CD 

II 

O 




Substitution of equations (A57 ) } (A 58 ), (A59 ) > and (A 60 ) into the 
homogeneous portion of the energy equation (A17) yields 


s A ax 1 dY _ _1_ ds\ = Cv^, d 2 0/dg 2 
X 3 Y 3 \X dx + Y dx " 2s dx/ f de/d^ 


(A62) 


If t, and x are independent variables, both sides of equation (a 62) must 
equal a constant, say -Cvoo 


s (± dX 1 dY _ 1_ ds\ 

X^Y 3 \X dx Y dx ” 2s d x/ ~ V ” 


(A 63 ) 


d g e 

d^ 2 


+ £“ f - 0 


(A64) 
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Since X is an arbitrary function of x, let 


1 dX l.dB =0 
X dx " 2 s dx 


(A65) 


Upon integration. 


X = C 3 s 1/2 


(A 66) 


Substituting equation (A66) into (A63 ) , there results 

fS = - Gv » C 3 3 s 1/2 (a 67) 

Upon integration with the boundary condition, x = xj, Y- 3 *- <», the arbitrary 
function Y is determined as 


Y = 



(a68) 


The variable t; is from equations (A60), (a 66), and (A68), 


£ = XYy 


s 1/g y 



(A69) 


The function S is found by integrating equation (A 64 ) with the boundary 
condition 0( 0) = 0. 


9 



(A 70 ) 


The constant C 4 is determined from the boundary condition, 0 (00) = 1. 
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C 4 



(A7U 


The complete solution of the homogeneous equation for the boundary 
conditions (A6l) is 


fl(S) 


-st — i— rK< 3 /^ 

%o-t s 3 i/3 r(^ 4 > 


(A72) 


General solution of energy equation .- Addition of solutions (A55) 
and (A72) with x j = x Q yields the temperature distribution in the 
boundary layer of a flat plate with the porous surface at a constant 
temperature T c and the nonporous surface at a different constant tem- 
perature T s . In order to obtain the solution when the temperature of 
the nonporous surface is arbitrary, homogeneous solutions of the form 
(A72) can be added as shown in reference 6, because the energy equation 
is linear in temperature . For r surface temperature discontinuities 
in the nonporous region, x > Xq, the temperature is given by 

j=r-i 

T(x,y) = T r T c+ £ (T Sj “T 8j+i )e(x,y, Xj ) + T Br (A73) 

j=l 


where 


t B;l = *c 


As r-> oo, the summation approaches an integral such that 


T(x,y) 


T I- T c +T s 


W 7 


dTsCx!) 

dxi 




(A7^) 


The general solution of the energy equation for an arbitrary surface 
temperature distribution in the nonporous region is then 


T(x,y) = 


T a (x)+(T t -T ( 


.) — - 

Uco 


u£ 


x 

- J e(x,y,x 1 )dx 1 (A75) 

Xq 
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It should be noted from equations (A9) and (A10) that 

dl s _ dT s » 

dx = dx* ~ (A 76 ) 


Seat transfer at the surface .- The heat transfer 
unit area is by definition and from equation (AlO) 


at the surface per 






Substitution of equation (A75) into (A77) yields 


(A77) 


<ls = 



x 0 



(A78) 


From equation (A 69 ) 



(A79) 


From equation (A72) 



(A80) 



Substituting equations (A79) and (A 80 ) into (A 78 ) 
for the heat transfer at the surface for a given 
distribution. 


results in the equation 
surface temperature 
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The variables are normalized by letting 


I 






(A82) 


Also, the following substitutions are made 


(A03) 

e 

<P(S) = / 'feti) dg 

Then, the heat transfer at the surface, equation (AfiL), in terms of the 
normalized variables becomes 


q«(e) - - •^ S2££ ® 




(Tt-Tc)^O) -0.53&> o 1/3 f 


6 


a *JL 


C«p(6)-q>(Si)3 


X/3 


dg: 


(A84) 


Surface temperature with prescribed heat transfer at the surf ace . - 
With the heat transfer at the surface prescribed, the surface temperature 
is given by the following integral equation based on equation (A 8^) . 


/ 


| dTs(ji) 

Mjl 

[<p(l ) _c p(Sx)3 1/3 


dg x = o(g) 


(A 85 ) 


where 


Bo Lfbo^CpN^I) 


(T t -T c )*s/^(g) 


(A86) 
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f r T hp. independent variable is now changed from 5 to 9 in equation (A 85) 
to yield 

dJs ( ) 


I 


<P 


(43= 


(AflT) 


The integral equation now has the same form as Volterra’s equation of the 
first kind with the kernel becoming infinite at the upper limit. Solution 
of this equation is given in reference 7 • 


dT s (cp) sin (jt/3) & 


dcp 


dcp 


/ 


9 0(<Pi)4«Pi 
(<P-9 i ) £/3 


(A88) 


Upon integration frcm cp = 0 to cp. 


T s (cp) = 


sin (it/ 3) f 9 fl(<Pi)d.<Pi 
* l (CP-93.) 273 


(A89) 


Now the independent variable is changed back to | 


_ sin (n/3) P fl(ll)N/p(gl) d6x 

sU) c * J [cp(D-cp(a)] E ' 3 


(A90) 


Consider the special case where the nonporous surface is insulated or 
qg(|) = 0 in equation (A86). Then 


0(6) 



(T t - T c)^(iT 


The temperature of the insulated surface is then 


T s (|)-T c = 0.512 f g(6i)d6i 

T t -T c D 0 1/3 J [9(6)-9(ei)] S,S 


(A91) 


(A92) 


Solution of equation (A92) was performed on an IBM 650 computing machine. 
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TABLE I.- EFFECT OF UPSTREAM TRANSPIRATION COOLING RATES 
ON INITIAL CONDITIONS FOR PRESENT ANALYSIS 


F=f (0) 

Do 

cr(l) 

0 

31.76 

1.871 

-.25 

37-77 

1.504 

-.50 

45-97 

1.115 

-•75 

58.83 

.718 

-1.00 

83.87 

.325 











Figure 4.- Comparison of seventh-degree polynomial velocity profile with Low's exact solution 

for various transpiration cooling rates. 
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Figure 8.- The effect on the temperature distribution of varying the 
porous length of a partially transpiration-cooled plate subject to 
constant total coolant flow at a fixed temperature . 
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Figure 8.- Concluded. 
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